function nodes = getnodes(m,lb,ub)
% chebyshev nodes
z    = -cos(3.14159*(2*(1:m)'-1)/(2*m));
nodes  = (z+1)*(ub-lb)/2+lb;
